A Nonperturbative Perspective on Inner Product Quantization: Highly Accurate 

Solutions to the Schrodinger Equation 

C. J. Tymczak, G. S. Japaridze, C. R. Handy and Xiao-Qian Wang 
Department of Physics & Center for Theoretical Studies of Physical Systems, Clark Atlanta University, Atlanta, Georgia 30314 

(Received February 9, 2008) 



We devise a new and highly accurate quantization pro- 
cedure for the inner product representation, both in config- 
uration and momentum space. Utihzing the representation 
= ^ . ai[i?]C'^/3(C)i for appropriate reference func- 
tion, RpiS,), we demonstrate that the (convergent) zeroes of 
the coefficient functions, ai[E] = 0, approximate the exact 
bound/resonance state energies with increasing accuracy as 
i — » oo. The vahdity of the approach is shown to be based on 
an extension of the HiU determinant quantization procedure. 
Our method has been apphed, with remarkable success, to 
various quantum mechanical problems. 

PACS numbers: 02.30. Hq, 03.65.-w, 03.65.Ge 

The development of numerically efficient algorithms for 
quantum mechanical problems continues to be a challeng- 
ing issue. Within this context, variational approaches 
based on optimal selection of basis sets remain one of 
the more successful approaches. Nevertheless, it appears 
that little attention has been paid to further enhancing 
such techniques by utilizing basis functions with an in- 
herently analytical (i.e., power series) significance. In 
this Letter we show that such considerations can greatly 
reduce the computational effort of variational procedures 
by significantly reducing the size of the problem. Many 
important quantum mechanical problems, heretofore in- 
volving large variational matrices, reduce to solving for 
the roots of an algebraically generated function. Our ap- 
proach is based on a relatively simple extension of the 
Hill determinant approach, which nevertheless, has not 
been exploited until now. 

One of the basic procedures for solving quantum sys- 
tems is the Hill determinant approach [1]. It involves 
approximating the bound state wave function in terms 
of a suitable truncated basis, 'i'{x) = "^I^q ViBi{x) , and 
solving the finite dimensional problem 

/ 

Y,M.,AE]vj^O, (1) 

through the Hill determinant equation 

Det(^M^^^[E]^ ^0, (2) 

= M^,, for < I, J < I, SindM^,[E] = {B,\H\Bj)~ 
E{Bi\Bj). For a suitable basis, as / — > oo, the roots of 
the Hill determinant converge to the true eigenvalues of 
the Hamiltonian. 



Another basic analytical tool is the power series expan- 
sion method, specifically in the form of the inner product 
representation [2]: 

^ix) = {J2<E]x')Rp{x), (3) 

i 

for an appropriate reference function Rfj{x). The func- 
tion coefficients, ai[E], depend on the energy variable, E, 
and are readily generated through the standard methods 
for differential equations. For simplicity, the above ex- 
pansion assumes that x = is a regular point. If the 
functions a;*i?/3(x) define a complete basis (not necessar- 
ily orthonormal) , one can also pursue a Hill determinant 
analysis for the corresponding representation 

^'(a;) =^i;,a;*i?0(x). (4) 

i 

In some cases, the recursive structure of the Hill de- 
terminant for increasing values of I can be computed 
[3]. This allows one to analyze the asymptotic behavior 
with respect to /, for the roots of the Hill-determinant 
equation. In general, this analysis can be difficult and 
computationally demanding. It is in this context that we 
have discovered a remarkable relation whose simplicity 
has apparently gone unrecognized until now, despite the 
suggestive, but specialized, nature of the work by Bender 
and Dunne [4,5]. Specifically, we demonstrate that the 
convergent zeroes of the coefficient functions: 

a,+i[£;f'] =0, (5) 

(where I lables the roots) converge to the exact discrete 
state energies, ^]^g expansion order, i, in- 

creases: 

hm Ef^ = . (6) 

i — ^oo 

The theoretical justification for this proceeds as fol- 
lows. Assume that Bi{x) = x^R^lx) and that the corre- 
sponding Hill determinant method yields convergent re- 
sults to the physical energies and wave functions. Let E 
assume any value, E — Ec, for which the infinite matrix 
■Mij[Ec\, has no minor sub-matrix with zero determi- 
nant. One can recursively generate [6], through an effec- 
tive LU decomposition method, an infinite set of vectors 
{V^''>\0 <I<oo} satisfying, 



1 



j=0 



(7) 



or 



for < i < 7 — 1, and 

I 



j=0 



(8) 



where = 1, and vj'^^ = 0, for j > I +1. One also has 

Det(^M'-^^[Ec]j =Ul^gV^[Ec]. The relatfon in Eq. (7) 

involves / constraints for I unknowns (recall f/^'' = 1, 
thus Eq. (7) is actually an inhomogeneous relation). The 
second relation, Eq. (8), serves to define VjIEc]. 

For a given order /, the roots of Eq. (2) corre- 

dcfined 

)■ 



(I) 
I ' 



sponds to the roots of Eq. (8), E = E^ 
by Det(^M^'<'\El^^]^ ^ and Det{MS^\E\^ 
T)i\Ef\ = 0. We denote the corresponding vectors by 

Prom the recursion formula for the F's [6] , we have: 



Vr'\EP] = 



pin 



(9) 



Vi[eP] 

Thus, in the E EY' limit, for a given one obtains 

ya+i) = ±00 _ (10) 

By the very nature of the Hill determinant approach, 
for a sufficiently large order, /+!, one expects the partial 
sums 



(11) 



define good, approximate, solutions to the physical wave 
functions. The same should hold for the corresponding, 
renormalized, truncated power series expansion 



(12) 



Note that each of the expansions is normalized by setting 
the highest degree coefficient to unity. 

For fixed, and sufficiently large, values of / + 1, the 
preceding partial sums are expected to agree, term by 
term, for each of the roots E^^^^\ In particular, based 



on this concurrence, one expects °"^^fL 



which upon comparing with Eq.(lO) leads to the desired 
conclusion: 



(13) 



a,+i[Ef^]=0, 



(14) 



and <l>i{x)^NY:U(>ai[E\%'R0{x). 

We now demonstrate the capabilities of the preced- 
ing method. For the case of exactly solvable models, 
e.g., V{x) = or V{x) = + g/x'^, our approach 
readily yields the exact solutions, once the proper ref- 
erence functions are selected: Rp ~ exp(— a;^/2), and 
i?/3 = x" exp(-a;V2) {a = + + 4:g)/2), respectively. 

Let us now consider the quartic anharmonic oscillator, 
V{x) = x'^ + gx'^. Using i?^ = exp(— /3a;^) one obtains 
the recursion relation: 



UniE) = 



0„ an-2{E) + (1 - 4/32)a„_4(£;) + gan-6{E) 



n(n — 1) 



(15) 



where = 4/3n — &(3 — E and a„ = for n < 0. Table I 
shows the calculated energies of the ground and first ex- 
cited states. Our method shows systematic convergence 
for increasing /, exceeding some of the high accuracy so- 
lutions published [6 - 8] . Figure 1 shows the dependence 
of the ground state energy on the parameter g. Figure 2 
shows the corresponding wave functions. 

An important version of the quartic anharmonic os- 
cillator potential is the double well problem V{x) = 
— Z^x^ + x''. It is well-known that in the deep well limit 
(Z^ oo), the two lowest states are almost degenerate. 
Application of our method (refer to Table II) readily con- 
firms this, and by its high accuracy nature, significantly 
disagrees with the predictions of de Saavedra and Buen- 
dia (SB) [9]. In particular, for Z'^ = 25, we observe 
that the quasi-degenerate nature of the ground and first 
excited state energies becomes apparent only after 26 sig- 
nificant digits, not the 16 predicted by SB. 

Similar calculations have been carried out for higher 
degree potentials such as the sextic, octic, and dectic 
anharmonic potentials. The results are given in Table 
III. 

The generality of our method permits the study of 
large classes of problems. In particular, transcenden- 
tal potentials can be analyzed, provided the potential 
function, V{x), admits a power series expansion which 
is monotonically convergent (non- alternating). For in- 
stance, in the case of V{x) = exp(a;^) — 1, we can readily 
obtain the first three energy levels: Eq = 1.356371240, 
El = 4.633078503, and E2 = 8.9706782. {R^ix) = e"^", 
and a„ generated up to n < 80). 

Another type of potential which can be investigated is 
V{x) = Vo{xo + x)", where a > —2. For the two cases 
— 2 < a <0 and < a < 00 the asymptotic forms for the 
wave function can be studied via '^{x) = (p-ix) e~^/~^^ 
and 5'(x) = (p+ix) e^^^ , respectively. Upon making the 
transformation y = -^log(l +x/xo), one can readily ob- 
tain the corresponding differential equations for 4'±{y), 



2 



and proceed to generate the necessary power scries ex- 
pansion, (j)±{y) = J2^n^^y"'- Figure 3 we plot E[a] vs. 
a for the cases of {xq = 1, Vq = —2; — 2 < a < 0} and 
{xo = 1, Vo = 1; < a < 4}. 

The selection of the reference function is clearly im- 
portant. In the case of the sextic anharmonic oscillator 
{V{x) = -V gx^\ our method works for the reference 
functions e~^'^ , a = 2, and 3. For the case a = 4 
(and /3 = which corresponds to the asymptotic form 
of the wave function, no convergent roots were observed 
upon solving Eq. (14). 

Quantization is a global problem, not a local one. 
Configuration space analysis, by its manifestly local na- 
ture, may not always be a suitable representation in 
which to quantize. However, a momentum space rep- 
resentation is more appropriate. Specifically, the power 
series expansion for the Fourier transform of the wave 
function, ^(fc) ^ Jdk e-^'-^ix) = Ep^o^^^b), 
involves global quantities, the (Hamburger) power mo- 
ments, /i(p) = j dx xP^(a;). As such, one would expect 
that our method would be more effective in momentum 
space. These observations are at the heart of various mo- 
ment based quantization approaches [2] , particularly the 
Eigenvalue Moment Method (EMM) [10]. 

Extending our method into momentum space presents 
additional features not encountered in configuration 
space. The most important of these is that more vari- 
ables (the missing moments within the EMM approach) 
are encountered, regardless of the spatial dimension, V, 
of the problem. Solving spatial problems of dimension 
D = 1 presents similar theoretical/algorithmic challenges 
to those in higher dimensions. We outline the essentials. 

Limiting ourselves to symmetric ^(a;) configurations, 
for simplicity, the coefficients of the power scries expan- 
sion for the momentum space wave function, 4'(/c) = 

/ dk e~^''^'i{x) = Ep=o ^^2p)\' ''^(P') satisfy a linear, mo- 
ment recursion equation, resulting in the linear relations 
"(P) = ET=oME{p,e)u{£), < p < oo; the M£(p,£)'s 
are known, and the missing moment order, m^, is prob- 
lem dependent [10]. 

Implementing our quantization procedure on the rep- 
resentation \&(/c) = (En=o '*2n^^")e~'^'^ ) we obtain 

a2n[E,u{0),...,u{ms)]=Y,Dr,,e[EHe), (16) 



where 



DnAE] 



E 

Pl+P2=n 



(2pi)!p2! 



(17) 



In principle, there will be a sequence of energy and miss- 
ing moment values satisfying a2„[-E''"\ (i')}] = 
converging to the respective physical state values as 
n ^ 00. Since the D„,^[i?]'s are not expected to de- 
fine a degenerate matrix for all E's, as n — > oo, we can 



approximate the converging sequence by considering the 
[nis -|- 1] X [nis + 1] matrix equation 



^2=0 



0, 



and the ensuing determinant equation. 



L>ei (£»(")[£;]) =0. 



(18) 



(19) 



Implementation of this for the quartic {mg = 1) and 
sextic (mg = 2) anharmonic oscillators yielded results 
consistent with those cited in Tables I and III. 

Some problems can involve no missing moments. One 
of these is the aforementioned sextic anharmonic oscilla- 
tor {nis = formulation), provided one implements the 
above formalism with respects to the configuration space 
expression ^{x) = '9{x)e^^^ . 



The same holds for the problem V{x) 



provided ^{x) 



2^ 



[11]. Table IV summa- 



rizes our results for this case, which surpass the accuracy 
calculated through an analytic continuation quantization 
procedure [12]. 

In summary, we have developed a nonperturbative 
approach for the inner product quantization procedure 
which allows one to calculate the energies and wave func- 
tions of the Schrodinger equation. We have demonstrated 
that our method yields excellent numerical results for 
various quantum mechanical problems. A full account 
of the approach, including the application to higher di- 
mensional models and resonant states, will be published 
elsewhere. 
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TABLE L The calculated ground and first excited state 
energies for the quartic anharmonic oscillator with g = 1. 
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FIG 1. The calculated ground state energy for the 
quartic anharmonic oscillator. 

FIG 2. The calculated ground and first excited state 
wave functions for the quartic anharmonic oscillator. 

FIG 3. The calculated ground and first excited state 
energies for the potential V{x) = Vo{x + Xq)"' for —2 < 
a<2. 



TABLE H. The calculated ground and first excited state 
energies for the potential V(x) = —Z'^x^ + x'^. 



Z'^ Parity E± 






+ 


1.060 


362 


090 


484 


182 


899 


647 


046 


016 






3.799 


673 


029 


801 


394 


168 


783 


094 


188 


1 


+ 


0.657 


653 


005 


180 


715 


123 


059 


021 


723 






2.834 


536 


202 


119 


304 


214 


654 


676 


208 
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+ 


-3.410 


142 


761 


239 


829 


475 


297 


709 


653 






-3.250 


675 


362 


289 


235 


980 


228 


513 


775 


10 
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-20.633 


576 


702 


947 


799 


149 


958 


554 


634 






-20.633 


546 


884 


404 


911 


079 


343 


874 


899 


15 


+ 


-50.841 


387 


284 


381 


954 


366 


250 


996 


515 






-50.841 


387 


284 


187 


005 


154 


710 


149 


735 


25 


+ 


-149.219 


456 


142 


190 


888 


029 


163 


966 


538 






-149.219 


456 


142 


190 


888 


029 


163 


958 


974 



TABLE HL The calculated ground state energies of the 
sextic, octet and dectic anharmonic potentials for g = 1. 



\-{,r) 


i.0 (l?ef. 


£() 


x' + x" 


1.435 624 619 


1.435 624 619 003 392 231 569 


x^ + x^ 


1.491 019 895 


1.491 019 895 662 






1.546 263 512 6 



TABLE IV. The first four symmetric state energies for the 
rational fraction potential V{x) = x"^ + ^_f^^; for A = g = 0.1 
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Figure 2, Submitted to Phys . Rev. Lett. 
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Figure 3, Submitted to Phys . Rev. Lett. 
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